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Abstract 



We study a large class of critical two-dimensional Ising models namely critical 
Z -invariant Ising models on periodic graphs, example of which are the classical 
Z^, triangular and honeycomb lattice at the critical temperature. Fisher |Fis66| 
introduced a correspondence between the Ising model and the dimer model on a 
decorated graph, thus setting dimer techniques as a powerful tool for understand- 
ing the Ising model. In this paper, we give a full description of the dimer model 
corresponding to the critical Z-invariant Ising model. We prove that the dimer 
characteristic polynomial is equal (up to a constant) to the critical Laplacian char- 
acteristic polynomial, and defines a Harnack curve of genus 0. We prove an explicit 
expression for the free energy, and for the Gibbs measure obtained as weak limit 
of Boltzmann measures. 



1 Introduction 

In |Fis66] . Fisher introduced a correspondence between the two-dimensional Ising model 
defined on a graph G, and the dimer model defined on a decorated version of this graph. 
Since then, dimer techniques have been a powerful tool for solving pertinent questions 
about the Ising model, see for example the paper of Kasteleyn [Kas67j . and the book 
of Mc Coy and Wu |MW73j . In this paper, we follow this approach to the Ising model. 

We consider a large class of critical Ising models, known as critical Z -invariant Ising 
models, introduced in |Bax89j , studied in |Mer01[ ICS06| . More precisely, we consider 
Ising models defined on graphs which have the property of having an isoradial embed- 
ding. We suppose that the Ising coupling constants naturally depend on the geometry 
of the embedded graph, and are such that the model is invariant under star-triangle 
transformations of the underlying graph, i.e. such that the Ising model is Z -invariant. 
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We suppose moreover that the couphng constants are critical by imposing a gener- 
alized self- duality property. The standard Ising model on the square, triangular and 
honeycomb lattice at the critical temperature are examples of critical 2^-invariant Ising 
models. 

Let G be an infinite Z^-periodic isoradial graph with critical coupling constants on the 
edges. Then, by Fisher the Ising model on G is in correspondence with the dimer model 
on a decorated graph Q, with a well chosen positive weight function v on the edges. The 
decorated graph Q is also infinite and Z'^-periodic. A natural exhaustion of the infinite 
graph Q by toroidal graphs is given by {Qn}n>i, where Gn = {V{Gn), E{Qn)) = Q/nZ'^. 
The dimer model on Qn is defined as follows. A dimer configuration M of Qn is a subset 
of edges of Qn such that every vertex is incident to exactly one edge of M. Let M.{Qn) 
be the set of dimer configurations of Qn- The probability of occurrence of a dimer 
configuration M is given by the dimer Boltzmann measure, denoted Vn'- 



Vn{M) 



HegAf 
2'n 



where Zn = YliM(^M{gn) IleeM ^e. is the dimer partition function. 

The first theorem of this paper proves an explicit expression for the free energy, and 
for the Gibbs measure obtained as weak limit of Boltzmann measures, of the dimer 
model on the infinite graph Q. It can loosely be stated as follows, a precise statement 
is given in Theorem [6] of Section HI Let K he, a, Kasteleyn matrix of the graph Q. It is 
a weighted, oriented adjacency matrix of the graph Q, defined in Section HTTl 

Theorem 1. 

1. The free energy of the dimer model on Q is: 

- lim \\og{Zn) = -TJT^—yj I f ^ogdet{K{z,w)) — — , 

where K{z,w) is the Fourier transform of the infinite Kasteleyn matrix K (it is 
a matrix of size V{Qi) x V{Qi)). 

2. The weak limit of the Boltzmann measures Vn defines a translation invariant 
ergodic Gibbs measure V . The probability of occurrence of a subset of edges ei = 
uivi, ■ ■ ■ ,em = UmVm of Q, in a dimer configuration of Q chosen with respect to 
the Gibbs measure V, is: 

neir • • > e-) = Pf ((i^-')ei,.. ,e™), 

where is the inverse of the infinite Kasteleyn matrix defined in and 

(i^~^)ei,... is the submatrix of K~^ , whose lines and columns are indexed by 
the vertices defining the edges ei, • • • ,em- 



2 



The statement of Theorem [T] is shnilar to Theorems 3.5 and 4.3 of [ KUS06j . see also 
[CKPOlj . obtained for dimer models defined on infinite Z^-periodic bipartite graphs. 
These results do not apply as such to our setting, where the graph Q on which our dimer 
model is defined, is not bipartite. Here is a description of the additional problems we 
had to solve when dealing with a non-bipartite graph. 

1. A key requirement in proving such a theorem is to have a precise description of 
the zeros of det{K{z, w)). In the bipartite case, Kenyon, Okounkov and Sheffield, 
prove that this zero set is a Harnack curve, which implies that there are at most 
two zeros on T^. Quite surprisingly, we also prove that in our case, the zero set 
is a Harnack curve, by relating it to the zero set of the Laplacian. This is the 
content of Theorem [2] below. 

2. In order to prove weak convergence of the Boltzmann measures, Kenyon, Ok- 
ounkov and Sheffield prove weak convergence of the measures on a sub-sequence 
of n's, which is quite easy once a description of the zero set of det{K{z,w))) is 
available, and use the convergence theorem for Gibbs measures of Sheffield [She05J 
to deduce convergence for every n. When the graph is not bipartite, Sheffield's 
theorem does not hold. We give a direct argument proving weak convergence of 
the Boltzmann measures for every n. Our techniques are general and robust, so 
that they could be used for dealing with dimer models on other non bipartite 
graphs. 

The description of the zero set of det{K{z,w)) is given by Theorem [2] below. An 
equivalent statement is given in Theorem [8] of Section [Sj where it is also proved. Let A 
be the Laplacian matrix of the infinite Z^-periodic isoradial graph G, and let A(z, w) 
be its Fourier transform. 

Theorem 2. There exists a constant c ^ 0, depending on the isoradial embedding of 
the graph G, such that: 

dei{k{z,w)) = cdet(A(z,u;)), 

and {{z,w) : det{K {z , w)) = 0} is a Harnack curve of genus 0. Its unique point 

on is (1, 1), and it has multiplicity 2. 

The main steps of the proof are: 

1. Show that the zero set of det{A{z,w)) is the spectral curve of the dimer model 
on an isoradial bipartite graph. The main tool to do this is an interpretation 
of detA(z,ti;) in terms of cycle-rooted spanning forests, which is an analog of 
Kirchhoff's matrix-tree theorem |Kir47j . Then use the result of Kenyon and 
Okounkov [KO06j by which spectral curves of dimer models on bipartite isoradial 
graphs are Harnack curves of genus 0. 
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2. Show that det{K{z,w)) = cdet{A{z,w)). First prove that det{K{z,w)) divides 
det(A(z, w)) by using an exphcit parametrization of the zero set of det(A(z, w)), 
then use the interpretation of det(A(z, if)) as cycle-rooted spanning forests to 
show inclusion of the Newton polygons. 

Outline of the paper 

Section 2: Definition of the critical Z-invariant Ising model. 

Section 3: Fisher's correspondence between the Ising and dimer models. 

Section 4: Study of the critical dimer model corresponding to a critical Z-invariant 
Ising model on a Z^-periodic graph. Definition of the infinite Kasteleyn matrix K. 
Proposition [5] proves the existence and uniqueness of the inverse of the matrix K. 
Statement and proof of Theorem [6l 

Section 5: Statement and proof of Theorem [51 

2 The critical Z-invariant Ising model 

Consider an unoriented finite graph G = {V{G), E{G)), together with a collection of 
positive real numbers J = {Je)e&E{G) indexed by the edges of G. The Ising model on 
G with coupling constants J, is defined as follows. 

A spin configuration o" of G is a function of the vertices of G with values in { — 1, +!}■ 
The probability of occurrence of a spin configuration a is given by the Ising Boltzmann 
measure, denoted : 



where 




is the Ising partition function. 

2.1 The Z-invariant Ising Model 

Z-invariant Ising models are Ising models defined on a class of embedded graphs which 
have the property of being isoradial. The corresponding coupling constants are natu- 
rally related to geometric features of these graphs, and satisfy the star-triangle relation 
defined below. 
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2.1.1 Isoradial graphs 



A graph G is said to be isoradial [Ken02t IKS05] . if it has an embedding in the plane 
such that every face is inscribed in a circle of radius 1, and all circumcenters of the 
faces are in the closure of the faces. Prom now on, when we speak of the graph G, we 
mean the graph together with a particular isoradial embedding in the plane. Examples 
of isoradial graphs are given in Figure [D below. 




Figure 1: Examples of isoradial graphs: the square lattice (left), the honeycomb lattice 
(center), and a more generic one (right). Every face is inscribed in a circle of radius 1, 
represented in dashed lines. 



To such a graph is naturally associated the diamond graph, denoted by , defined 
as follows. Vertices of G^ consist in the vertices of G, and the circumcenters of the 
faces of G. The circumcenter of each face is then joined to all vertices which are on the 
boundary of this face, see Figure[2j Since G is isoradial, all faces of are side-length-1 
rhombi. Moreover, each edge e of G is the diagonal of exactly one rhombus of G*; we 
let be the half angle of the rhombus at the vertex it has in common with e, see 
Figure El (right). 




Figure 2: An isoradial graph (left). The white dots are the circumcenters of the faces, 
that are also the vertices of the dual G* . Its diamond graph is represented in the center. 
On the right is the half-rhombus angle 6^ assigned to an edge e. 
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The same construction can be done for toroidal isoradial graphs in which case the 
embedding is on a torus. 



2.1.2 Star-triangle transformation 

In this section, we suppose that G is a toroidal isoradial graph. It is then natural to 
choose the coupling constants J of the Ising model on G, to depend on the geometry of 
the embedded graph: let us assume that J is a function of 9^, the rhombus half-angle 
assigned to the edge e. 

The star-triangle transformation is a local operation on the graph G which preserves 
isoradiality. If G has a vertex of degree 3 (star), it can be removed and replaced 
by a triangle reconnecting its neighbors, see Figure O The Ising model is said to be 
Z -invariant if the Boltzmann probability of any event not involving the spin of the 
degree-3 vertex is preserved under this operation. 



Figure 3: The star-triangle transformation corresponds to a rearrangement of three 
rhombi of the diamond graph. Isoradiality is preserved by this transformation. 

Imposing this condition is a strong constraint. Indeed, there is only a one parameter 
family of coupling constants which defines a Z-invariant Ising model: J{9e) must satisfy 
the equation: 



The functions sn, cn are Jacobi elliptic sine and cosinqj, and K is the complete integral 
of the first kind: 





(1) 




A detailed derivation of this parametrization is given in |Bax89j . Chap. 6. 



See [RG57] . Section 6.1 for definitions and properties of these functions. 
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2.2 Polygonal contour expansions of the Ising partition function 

Let G be a toroidal graph, with couphng constants J, not necessarily isoradial. There 
are two ways to define a probability measure on polygonal contours from the Ising model 
on G: the so-called low-temperature and high-temperature expansions of the partition 
function. 

The high temperature expansion is based on the following identity: 

f,J^TuCy ^ cosh( Je) (1 + tanh( Je)crucr^) . 

The partition function can be rewritten as 

Z^=[\{ COSh(Je)j ^ JJ (l + C7„(7„tanh(Je)) 

^e&E{G) ' ae{-l,l}^(^'' e=™eE(G) 




where V is the family of all polygonal contours drawn on G, for which every edge of G 
is used at most once. 

The low-temperature expansion gives a probability measure on polygonal contours on 
the dual graph G* of G, which represent the separations between the clusters of + spins 
and those of — spins. This expansion is based on the following identity: 

The partition function is thus equal to 




where V* is the family of polygonal contours drawn on G* , for which every edge of G* 
is used at most once. 

In both cases, the probability of a polygonal contour configuration is defined to be 

proportional to the product of the weights of the edges it contains. The weight of 
an edge e is tanh(Je) in the high temperature expansion, and exp(— 2Je.) in the low 
temperature expansion. 
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2.3 The critical Z-invariant Ising model 



Consider a Z-invariant Ising model on a toroidal isoradial graph G, and a Z-invariant 
Ising model on its dual graph G* . Let J and J* be the associated coupling constants 
satisfying ([1]) for some parameters k and k* respectively. Then, if k* = -7==, the 
measure on polygonal contours of G obtained from the high temperature expansion 
on G, is the same as the measure on polygonal contours of G obtained from the low 
temperature expansion on G*: there is a duality relation between Z-invariant Ising 
models with parameters k and k* . It differs from the classical one on the square 
lattice jKW4ij in the sense that the graph G and its dual are a priori not related: this 
duality relation must be viewed on the class of all Z-invariant Ising models. 

The only value of k for which the set of Z-invariant Ising models with parameter k 
maps onto itself (that is A;* = /c), is /c = 0. At this special value, we say that the model 
satisfies a generalized self-duality property. The expression for the weights simplify 
drastically, indeed: 



K{{)) = -, sn - 



2' \ TT 

The expression of the coupling constants reduces to: 



COS ( 



j{e,) = -\og[———]. (4) 



2 \ COS 6e 

In the case of the square lattice, the rhombus half-angle assigned to each edge is 7r/4, 
so that the coupling constants are all equal to log \/l + ^/2, which is the value of the 
critical temperature of the Ising model on the square lattice |KW41] . Similarly, on the 
honeycomb (resp. triangular) lattice, the rhombus half-angle assigned to each edge is 
7r/3 (resp. tt/6), in this well the coupling constants coincide with the critical 

temperature log \/2 + \/3 (resp. log \/ V3) [Wan45j . For these reasons, we refer to 
A; = as the critical value of the parameter, and the choice of weights ^ as the critical 
coupling constants. 

The Z-invariant Ising model on an isoradial graph with this particular choice of coupling 
constants is referred to as the critical Z-invariant Ising model. 



3 Fisher's correspondence between the Ising and Dimer 
models 

Let G be any graph drawn on a torus, with coupling constants J. Then, Fisher |Fis66| 
exhibits a correspondence between the Ising model on the graph G, and the dimer 
model on a "decorated" version of G. This correspondence uses the high temperature 
expansion of the Ising partition function given in Let us first recall the definition 
of the dimer model. 
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3.1 Dimer model 



Consider a finite graph G = {V{G), E{G)), and suppose that edges of Q are assigned a 
positive weight function u = {i^e)eeE{g)- The dimer model on Q with weight function v 
is defined as follows. 

A dimer configuration M ofQ, also called perfect matching, is a subset of edges of G such 
that every vertex is incident to exactly one edge of M. Let M{G) be the set of dimer 
configurations of the graph Q. The probability of occurrence of a dimer configuration 
M is given by the dimer Boltzmann measure, denoted V^: 



where Z'^ = J^MeMiG) IleeAf is the dimer partition function. 



3.2 Fisher's correspondence 



Consider an Ising model on a toroidal graph G, with coupling constants J. In this 
paper we use the following slight variation of Fisher's correspondence. 

The decorated graph, on which the dimer configurations live, is constructed from G 
as follows. Every vertex of degree of G is replaced by a decoration consisting of 3k 
vertices: a triangle is attached to every edge incident to this vertex, and these triangles 
are linked by edges in a circular way, see Figure H] below. This new graph, denoted by 
G, is also embedded on the torus and has vertices of degree 3. It is referred to as the 
Fisher graph of G. 




Figure 4: A vertex of G (left) with its incoming edges, and the corresponding decoration 
in Q (right). 

Here comes the correspondence: to any contour configuration C coming from the high- 
temperature expansion of the Ising model on G, we associate 2l^('^)l dimer configu- 
rations on G: edges present (resp. absent) in C are absent (resp. present) in the 
corresponding dimer configuration of Q. Once the state of these edges is fixed, there 
is, for every decorated vertex, exactly two ways to complete the configuration into a 
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dimer configuration. Figure [5] below gives an example in the case where G is the square 
lattice I?. 



1 2 



3 4 




Figure 5: Polygonal contour of 1? , and corresponding dimer configurations of the 
associated Fisher graph. 



Let us assign, to an edge e of ^, weight fe = 1, if it belongs to a decoration; and weight 
fe = coth(Je), if it corresponds to an edge of G. Then the correspondence is measure- 
preserving: every contour configuration C has the same number (2l^*^'^)l) of images by 
this correspondence, and the product of the weights of the edges in C, ngg^tanh(Je) 
is proportional to the weight He^c coth( Jg) of any of its corresponding dimer configu- 
rations with a proportionality factor, nee£;(G) tanh(Je), which is independent of C. 

Note that Fisher's correspondence is in fact valid for any graph G drawn on a surface 
without boundary, in particular, for an infinite planar graph. 



3.3 Critical dimer model on Fisher graphs 

Consider a critical Z-invariant Ising model on a toroidal isoradial graph G. Then, the 
dimer weights of the corresponding dimer model on the Fisher graph Q are: 

1 if e belongs to a decoration, 

coth(J(6'e)) = coth {log ^cos e!" ) = (%) ^ comes from an edge of G. 

We refer to these weights as critical dimer weights, and to the corresponding dimer 
model as critical dimer model on the Fisher graph Q. 

Because of the self-duality property of the critical Z-invariant Ising model, the same 
dimer model is obtained when considering Fisher's correspondence for the low temper- 
ature polygonal contour expansion of the Ising model on the dual graph G* . 

Remark 3. Critical dimer models on Fisher graphs do not enter in the framework of 
critical dimer models introduced in [Ken02] . Indeed, the latter are restricted to live 
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on bipartite isoradial graphs, whereas Fisher graphs are neither bipartite nor isoradial 
(although they are constructed from isoradial graphs). 

In Section [3] below, we consider critical dimer models on infinite Z^-periodic Fisher 
graphs. The main results are Proposition [5] proving existence and uniqueness of the 
inverse Kasteleyn matrix, and Theorem [6] giving an explicit expression for the free 
energy and for the Gibbs measure obtained as weak limit of Boltzmann measures. 



4 Critical dimer model on infinite periodic Fisher graphs 

Let Q = [V{Q), E{Q)) be an infinite Z^-periodic Fisher graph obtained from an infinite 
isoradial Z^-periodic graph G. Assume that edges of Q are assigned the dimer critical 
weight function denoted by i^. 

Consider the exhaustion {Qn}n>i of Q by toroidal graphs, where Gn = Q jnTl? . The 
graph Q\ is called the fundamental domain of Q. Let us assume that Gi has an even 
number of vertices, (if this is not the case, the graph Gi has no dimer configuration). 

In Section 14. 1^ we define Kasteleyn matrices and state Theorem H] which is the yet 
classical result proving an explicit expression for the partition function of Qn, denoted 
by Zn, and for the the Boltzmann measure of Qn, denoted by Vn- Theorem H] is valid 
in greater generality, i.e. for dimer models on Z^-periodic graphs with positive weight 
functions. We nevertheless choose to state them in the context which is relevant to this 
paper. 



4.1 Dimer model on toroidal graphs Qn 

The key objects used to obtain explicit expressions for the dimer model are Kasteleyn 
matrices. They are weighted oriented adjacency matrices of the graph Gn defined as 
follows. 

A Kasteleyn orientation of the edges of Gn is an orientation of the edges such that all 
elementary cycles are clockwise odd, i.e. when traveling clockwise around the edges of 
any elementary cycle of Gn, the number of co-oriented edges is odd. If the number 
of vertices of the graph is even, then such an orientation always exists |Kas671 [CR07) . 
The Kasteleyn matrix corresponding to such an orientation is a |T^(^„)| x |T/(^„)| skew- 
symmetric matrix, denoted Kn, defined by: 



u,v 



Uuv if n ~ f , and u ^ v 
—Vuv if n ~ f , and u ^ v 
else. 



Let 7a;^„ (resp. 7j/,n) be a path in the dual graph of Gn, winding once around the torus 
horizontally (resp. vertically). For 9,t & {0, 1}, let K^'^ be the Kasteleyn matrix K^, 
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in which the weights of the edges crossing the horizontal cycle 'yx,n are multiplied by 
(— 1)^, and those crossing the vertical cycle 7^^^ are multiplied by (—1)^. Then, it is 
always possible to choose the Kasteleyn orientation of the edges of Qn to be periodic, 
and such that TheoremfUbelow holds (changing the Kasteleyn orientation would merely 
change the order of the linear combinations below). 

Theorem 4. IKas61l OCPOTl lDZM+961 ITesOOl IUROTI . IKen97l 



1. The dinner partition function of the graph Qn is: 

Zn = \{- Pf «) + Pf(^^°) + Pf + ■ (5) 

2. Let E = {ei — "Uifi, ■ ■ ■ , Cm — UfYiVm} ; be a subset of edges of Qnj then the 
probability Vn{ei,- " -.^m) of these edges occuTvifig in. a diTficv coTifiguTcitiofi of Qn 
chosen with respect to the Boltzmann measure Vn, is: 



Pf(i^r + Pf(^r + PiiKlE^ + Pf(^„ , (6) 



where E'^ = V{Qn) \ {ui-,vi,--- ,Um,Vm}, and {K^)e'^ is the submatrix of Kf^ 
whose lines and columns are indexed by E^. 



4.2 Critical dimer model on infinite periodic Fisher graphs 

In this section, we start by defining the Kasteleyn operator of the infinite Z^-periodic 
Fisher graph then Proposition [5] shows existence and uniqueness of its inverse. The- 
orem [6] proves an explicit expression for the free energy, and for the Gibbs measure 
obtained as weak limit of Boltzmann measures for the critical dimer model on Q. 



4.2.1 Kasteleyn operator and characteristic polynomial 

Assume that the infinite graph Q has the periodic Kasteleyn orientation induced by 
the Kasteleyn orientation of Qi. Let K be the corresponding infinite Kasteleyn matrix. 
Then K is an operator acting on 

The operator K acts by convolution when interpreted as follows. Since Q is Z^-periodic, 
vertices of G can be identified with V{Gi) x 7?, i.e. {v, x, y) is the vertex v in the (x, y) 
copy of the fundamental domain Qi. Coefficients of the infinite Kasteleyn matrix K 
are thus indexed as K^y^x^y-^^(^yi ^yiy A function / G C^(^) can be written as a complex 
vector- valued function f :!? ^ C^(^i), 

{f{x,y))v = f{v,x,y). 
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Since K is translation invariant, it acts as: 

{v',x',y') {v',x',y') 

= {{K*f){x,y)),, 

where K : I? ^ M{\V{Qi)\ x is the complex matrix-valued function defined by: 

K{x,y)vy = ^(^,,^:,y),{^;',o,o)• 
Introducing a new notation for the matrix-valued function defined above seemed arti- 
ficial to us, we nevertheless hope that this will not confuse the reader. 

Let ^ : T2 ^ M{\V{gi)\ x \V{gi)\) be the Fourier transform of K: 

K{z,w)= K{x,y)z''wy. 

(a;,j/)eZ2 

The dimer characteristic polynomial of the graph Q, denoted by P{z,w), is defined to 
be the determinant of K{z, w). The set 

{iz,w) G : P{z,w) = det{K{z,w)) = 0}, 

is the dimer spectral curve of the graph Q. 

Theorem [21 see also Theorem [8] of Section [5] states that the dimer spectral curve of the 
graph ^ is a Harnack curve of genus 0, and that the dimer characteristic polynomial 
of the graph Q has a unique zero on T^, located at (1, 1), with multiplicity two. The 
proof of Theorem [8] is postponed until Section \5\ These properties of the characteristic 
polynomial of G are used in a crucial way in the next two sections. 



4.2.2 Inverse of the Kasteleyn operator 

Let us define the space of rapidly decaying vector- valued functions S{1?), and the space 
of vector- valued functions decaying at infinity Co{l?): 

S{I?) = {f -.1? ^ C^(^i) : V(m,n) E N^, lim y)|| = 0}, 

\\{x,y)\\^co 

Co(z2) = {/ : ^ C^(^i) : Ve > 0,3iV G N, sup ||/(x,y)|| < e}. 

\\{^,y)\\>N 

Proposition [5] below states the existence and uniqueness of the inverse Kasteleyn oper- 
ator. The proof is inspired from the classical one for the continuous Laplacian on R", 
using Fourier transform. 

Proposition 5. For every function f £ S{1?), there exists a unique function u G 
Cq{I?), such that Ku = f. 
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Proof of existence. ByTheorem[2l for {z,'w) E T^\{(1, 1)}, P{z,w) = det{K(z,w)) ^ 0. 
As a consequence, the following matrix valued function is well defined on \ {(1, 1)}: 



K {z,w) = [K{z,w)] = — ^ 



P{z,w) 

where Coi{K{z,w)) is the cofactor matrix of the matrix K{z,w). Moreover, since 
K{1, 1) is skew-symmetric of even size and not invertible, the dimension of its kernel 
is even and non zero, thus at least equal to 2. Therefore Coi{K{l, 1)) is identically 0, 
and the singularity of K~^{z,w) at (1, 1) is of order 1 at most, and not of order 2 as 
one could have expected. This implies that is integrable on T^. Its inverse Fourier 
transform is thus well defined, let us denote it by U: 



U{x,y)=J^-\k~^){x,y) 



Co{{k{z,w)y _^ _.^dzdw 



{2'KiY J Jf2 P{z,w) z w 

Define the complex vector- valued function u '.W? ^ (^yiGi) |-,y. 

u = U * f. 

Let us first prove that Ku = K * u = f . For every {x, y), (x2, 2/2) £ we have: 
^ K{x - xi,y - yi)U{xi - X2,yi - y2) = 

(a;i,s/i)eZ2 

Since the sum over is in fact finite, we deduce: 



{K *U * f){x,y) = ^ ^ K{x-xi,y -yi)U{xi-X2,yi-y2)f{x2,y2) 
(xi,j/i)ez2 (x2,2/2)ez2 

= ^ 6x^xjy^y^ld\v{gi)\ f{x2,y2) = f{x,y). 

Let us now prove that u is in Cq{I?). Split the sum over {x' ,y') into two parts, depending 
on the distance to {x,y): 

u{x,y) = ^ U{x-x ,y-y)f{x ,y)+ ^ U{x-x\y-y)f{x,y). 

{x',y') (x',y') 
d{{x,y),{x',y'))>L d{{x,y),{x' ,y'))<L 

By classical estimate on the Fourier coefficients of a rational fraction (see for example 
|KQS06| . Lemma 4.4), we know that U{x — x' ,y — y') decays linearly with the distance 
d{{x,y), {x' ,y')). We can therefore choose L large enough, so that 

\\U{x -x',y- y')\\ < 



+ r 
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when d{{x, y), {x' , y')) > L. The first sum is thus bounded by: 

^ \\U{x-x',y-y')\\\\f{x',y')\\< ' J] \\f{x',y')\\<'-. 
{x',y') {x',y') 

di{x,y),ix',y'))>L 

Since the coefficients U{x — x',y — y') are bounded and / G S{1?), one can choose 
= N{L) large enough so that if ||(x,y)|| > N, then for all {x',y') such that 
d{{x,y), {x' ,y')) < L, we have < — , where ul is the number of points 

in the ball of radius L. In that case, the second sum is thus bounded by: 



^ ||f/(x-x',y-yOllll/(^',y')ll <nL||t/||oo sup \\f{x',y')\\ < -. 

A{(xX(^'!y'))<L di{x,yUx',y')}<L 



Thus, when ||(x,y)|| > N, \\u{x,y)\\ < e, and we conclude that u S Co(Z^). 

Proof of uniqueness. Define ^(T^) to be the space of smooth functions on the torus with 
values in C^^^^^. Recall that the Fourier transform ^ is a bijection between S{7?) and 
5(T^). Let us denote by ( , ) the duality bracket between S{7?) and its dual S'{I?), 
and by ( , ) the one between 5(T^) and its dual 5'(T^). 

By duality, the Fourier transform extends bijection from ^'(Z^) to ^'(T^): 

Suppose that we have two functions ui and U2 in Co{Z,'^), such that K*ui = K*U2 = f ■ 
Let u = U2 — u\. The function u G Co{'I?), defines a distribution u in S'{1?). 

Then, since * n = 0, we have for every ip G 5(T^), 

= {J^{K*u),ij) = {K*u,J^-'^^) = -{u,K 

= -{u,J^-\Kij)) = -{J^{u),k^). (7) 

The third equality results from the fact that if is a skew-symmetric operator. 

Let be a function in 5(T^) whose support is contained in T'^ \ {(1, 1)}. Since K~^, 
as a function of {z,w) is smooth on \ {(1, 1)}, we deduce that 

^ = k~^^ G 5(T2). 

Plugging this in ([7]), we deduce: 

G cS(t2) with supp(e) c t2 \ {(1, 1)}, = 0. 

This implies that the support oi J^{u) is contained in {(1, 1)}. As a consequence, J^{u) 
is a linear combination of derivatives of the Dirac distribution at (1,1). But the only, 
such linear combination which has an inverse Fourier transform in Cq{1?) is the trivial 
one identically equal to 0, therefore, !F{u) = 0, i.e. u = 0. □ 
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In the light of Proposition \5\ let us define the inverse Kasteleyn matrix, denoted by 
K~^, to be the infinite matrix whose coefficients are: 

1 I ff Cof(j?(z,w)V,^ _^^,_^^ y-ydz dw 



4.2.3 Free energy and Gibbs measure 

The free energy of the dimer model on the infinite Z^-periodic graph Q is denoted /, 
and is defined by: 

/ = - lim ^log2:„. 

A Gibbs measure on the set of dimer configurations of ^, is a probability measure 

on Ai{G), which satisfies the following. If one fixes a perfect matching in an annular 
region of Q, then perfect matchings inside and outside of this annulus are independent. 
Moreover, the probability of occurrence of an interior matching is proportional to the 
product of its edges weights. 

In order to state Theorem [6] giving an explicit expression for the free energy of the 
dimer model on Q, and for the Gibbs measure obtained as weak limit of the Boltzmann 
measures, we need the following notation: ^ is the cr-field generated by cylinders, a 
cylinder being the set of dimer configurations of Q containing a fixed finite subset of 
edges of G. 

Theorem 6. 

1. The free energy of the dimer model on Q is: 

f = ~^T7^r~\2 [ [ ^ogP{z,w) — — . 
2{2-KiY J Jj2 z w 

2. There is a unique probability measure V on such that the probability 
of occurrence of a subset of edges E = {ei = uivi, ■ ■ ■ , = UmVm}, of G in a 
dimer configuration of G is: 

nei,--- ,em) = (n^"^.-^) Pf((^~'k), (9) 

where {K~^)e is the submatrix of the infinite matrix given in ([H]), whose 
lines and columns are indexed by the vertices defining the edges of E. Moreover, 
V is a translation invariant ergodic Gibbs measure. 

Proof. By Theorem [51 we know that P{z,w) has a unique double zero at (1, 1) on the 
unit torus. With this information, the general argument of [CKPOl] works for Part 1. 
We nevertheless sketch the proof for the convenience of the reader. 

For the proof of Parts 1 and 2 of Theorem [6], we need the following facts. 
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(a) This argument can be found in jCKPOlj . For every 9,t £ {0,1}, using Fourier 
transform, the matrix K^'^ can be block diagonalized, with (n^)^ blocks of size 

1^(^1)1 X \vig,)\: 

= {Qin^Bi^Qi^ (10) 

where the matrices Q^^ and are given in block form (x, y, j, /c, j'. A;' G { 0, • • • , n- 
1}): 

{Qnli;j,k),i;x,y) = -e' " ""e* n yid\v{g,)\, 
[Bn )i-,j,k),i-,j',k') = Oj,j'Ok,k'^ [e n ,e " j. 

Since the matrix B^"^ is block diagonal, we simply write for the (•, j, /c) 

diagonal block. The matrix Q^^ is unitary, i.e. {Q^)^Qfi^ = Id(|y(gj)|„2). 

(b) Recall that P{z,w) = det{K (z , w)) , and that by Theorem[2l P{z,w) has a single 
double zero at (1, 1) on the unit torus. We will use the following facts: 

o (e'—^,e'—^j = (1,1) ^ {j,k) = (0,0) and (0,r) = (0,0). 

/ . (2j + 9)7r . (2fc + T)7r \ 

o For all other indices, the points I e* « , e' « I are at distance at least 
O (^) from the zero (1, 1) of P on the torus. 

Proof of Part 1. 

As a consequence of Equation (fTOll : 

n— 1 n— 1 



det = det = n n ^ (-^"^ ' -''^) ■ (11) 

j=0 fc=0 

Using (b), we deduce that for every n, det(iC^'^) = 0, if and only if {0,t) = (0,0). 

Moreover, for {6, r) ^ (0, 0), the points (e* « , e* " ) are at distance at least O (-) 
from the zero of P, so that by the general argument of [CKPOlj (see also |KOS06| ). we 
have: 

lim \\ogV{{K'j) = I lim ^logdet«^) 



--, — ^ lim 


(2^) 


2 (27r)2 n-.oo 




1 1 /"^TT 




2(27r)2yo 





9 rt— 1 n— 1 

• .(2j + e)^ (2fc + r)7r 

e " , e 



j=0 fc=0 

logP(e^^e*^)d??c^C. 
By Theorem m 

max{Pf(K^)} < Z„ < 2max{Pf(iCf )}, (12) 
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and we deduce Part 1 of Theorem [6l 
Proof of Part 2. 

Let us prove that the Boltzmann probabihty Vni^ir " ^^m) given in Q of the edges 
ei = uivi,--- ,em = UjnVm being in a dimer configuration of Qn (n large enough) 
converges to the RHS of ([9j), as n ^ oo. The existence of a unique measure V equal to 
the RHS of on cylinder sets is then given by Kolomogorov's extension theorem (a 
similar argument is done in details in [dT07] ) . The Gibbs property results from the fact 
that P is a weak limit of Boltzmann measures. Translation invariance is a consequence 
of the periodicity of the weights, and ergodicity (even mixing) readily follows from the 
decay of the coefficients of with the distance between vertices. 

Case 1: Terms of ([6]) involving 'Pi{K^)E'^, {0,t) 7^ (0,0). In this case, the general 
argument of [ CKPOlj works (see also [KOS06] ). It runs as follows. By (b), for every n, 
det(-ftr^'^) 7^ 0. We thus apply Jacobi's formula for the determinant of minor matrices. 
Observing that all matrices involved are skew symmetric, this yields: 

pf(Kr)^c=pf((i^f)-i)^pf«-). 

Moreover, by (a), we have for all {6,t) ^ (0,0): 

^ n >{v,x,y),{v' ,x' ,y') 

^ n-ln-1 Cof K e " ,e 

j=o k=o -Tie " ,e 

/ ■ (2j + fl)7r . (2fc + T)7r \ 

Since the points I « , e* " I are at distance at least O (-) from the zero of P, 



the contribution of the points close to the singularity is negligible, and the Riemann 
sum converges to K'^^^y^^^^,^^,^^,^ given in ([H]). 

Moreover, ^^^p- is bounded, and since Pf(-fC°°) = 0, E(e,r)^(o,o) ^^2zP = 1- We 
conclude that: 

where {K~^)e is the sub matrix of the infinite matrix K^^, whose lines and columns 
are indexed by vertices of E. 

Case 2: Term of ([6]) involving Pf(K°'^)£;c. Let us show that, lim„^oo ^^^2^ = 0. 
We actually prove the following precise estimate: 

To simplify notations, we write = Kn, Q^^ = Q, B^^ = B. 



18 



By (a), B is skew-Hermitian, so that the diagonal block -B(.^o,o) ^oo- Moreover 
by (b), -B(.^o,o) = is not invertible since det(^(l,l)) = P(l, 1) = 0. Let 2q 

be the dimension of the kernel of -B(. 0,0)) ^i^d Pq be an orthogonal matrix of size 
X 11/(^1)1, such that: 

^(■,0,0) = -fo-^o-Po, 

where Dq is diagonal, with the first 2q elements being 0. Observe that Pq and Dq are 
independent of n. 

Define Q to be the matrix obtained from Q by multiplying all blocks of the first block 
line by P on the left. Then, Q is unitary and 

Kn = Q^BQ, 

where B is block diagonal, with: 

B(.,o,o) = -^0, 

^(■,j,k) = -B(.j-fc), for all {j,k) / (0,0). 

Let B be the matrix obtained from B by removing the first 2q lines and first 2q columns. 
Since, for {j,k) 7^ (0,0), all blocks B(^.j j^^ are invertible, the matrix B is. 

The proof of (fT3]) is a consequence of: 
Lemma 7. We have the following estimates: 



VdetB ^ V^^^ 

^ VdetB ^. , 
2. — = 0{n). 

Before proving Lemma [71 let us end the proof of Theorem [6l Combining the above two 
estimates, we see that the LHS of (fTSl) is 0{n^^^'^), and since g > 1, it goes to zero at 
least as fast as n~^, as n goes to infinity. □ 



Proof of Part 1 of Lemma^ Recall that \E\ = 2m. To simplify notations, let us set 
2N := |l/(^i)|n^. We will use the following convention: if A is a matrix, and S (resp. 
T) is a subset of row-indices (resp. column-indices) of A, the is the submatrix 
obtained from A by keeping rows with indices in S and columns with indices in T. If 
S (resp. r) is omitted, then all the rows (resp. columns) are keplH. 

■^To comply with the notations used before and to simply notations, there is one exception to this 
rule: {Kn)E'' still denotes the submatrix of K„ with rows and columns indexed by E"^ . 
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We apply Cauchy-Binet's formula to compute det{Kn)E'=' 

detiKn)E^ = det((QtBQ)|:) = det(Q^,,BQ^^) 

J2 det(Q^!)det(B^)det(Qf ). 

S,T C{1,- ,2N} 
\S\=\T\=2N-2m 

Since the first 2q lines and 2q columns of B are 0, det(ii'„) = if m < g, so that ([13]) 
is true. Let us assume that m > q. Then, 

det{Kn)E^= Yl det(Ql|)det(B^)det(Qf ). 

S,Tc{2q+l,- ,2N} 
\S\ = \T\=2N-2m 

Let S"" = {2q + , 2N} \ S, = {2q + 1, ■ ■ ■ , 2N} \ T. Then, by Jabobi's formula 

and unitarity of Q, we have: 



det(B^) = det(B)det((B-i)|c), 
det(Qf = det(Q) det(Qy^'), det(Q^"!) = det(Ql') det(Qfc 

So that: 



det((i^„)Bc) = det(B) J2 det(Q^'^') det((B-^)|c) det(Qfc). 

S,Tc{2q+l,- ,2N} 
\S\=\T\=2N-2m 

We have the following expansions for the determinants, see for example [Bou70| . A III.98, 
formula (21): 

det(Qfc)= PL,L'det(Qfi ... 2^|)det(Q|l), 

LCE 
\L\=2q 

det(Q^,^^)= Y PK,i^'det(Qj,^^'-'''^^)det(Qtf ), 

\K\=2q 

where L' = \ L, K' = \ K, and Pl,l'^ Pk,k' are signature of permutations 
constructed from L, L' and K, K' , see |Bou70j for a precise definition. 

Thus, 



det(i^„)£;c ^ A 4-/r\H^r-,2q}-, , ./r^L 

^g^(g) = 2^ PK,K' PL,L' det{Q'j^ ')det(Q|i ... 2,|) 

\K\ = \L\=2q 



Y det(QtT ) det((B-i)f:) det(Qg). 

S=T'=C{2g+l,-,2Ar} 
\S''\ = \T''\=2m-2q 
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Let Q be the matrix of size (2A^ — 2q) x 2N obtained form Q by removing the first 2q 
hnes. Then, using Cauchy-Binet's formula again, we have: 



det(Er, 



PK,K' PL,L' det(Qy^'- •'^^) det(Qfi,... det((QtB-iQ)|;.',)- 



det(B) 



K,LCE 
\K\ = \L\=2g 



Now, by definition of the matrix Q, all elements of Q|^^' '^''^ and Qi"-!^ 0{n ^). 

Hence, by Hadamard's inequality, det(Qy-^' '^'^^) and det(Q|'-^ 2q}) ^^'^ O {j^) ■ More- 
over, by definition of the matrix Q and B, 



(Q^B-iQ)(, 



{v,x,y){v' ,x' ,y') 



n-1 n-1 Cof LfC e " , e n 

j=o fc=o ( e " , e " 

(i,fc)7^(o,o) 



for some number c^^/ independent of n, coming from the contribution of the first block 
of B. Since {j,k) ^ (0,0), all points ^e*^,e*~^ are at distance at least O (^) from 

the of P, and the Riemann sum converges to K^J^ ^, given in dH]). Moreover, 
these terms are uniformly bounded, so that: 



det(QtB-iQ)f;, 



for some constant C. Since k and q are independent of n, we deduce that: 



det((ir„)j;c) 
det(B) 



□ 



Proof of Part 2 of Lemma^ Since the partition function Zn is bounded from below 
by any of the Pf(K^'^), it suffices to show that 



det(B) = det(i^^;^°) x O(n^). 



By (dH) we have: 



idet(oi= n 

i,fce{-LtJ,...,LtJ-i} 

|det(B)|=Ao n 

i,fce{-LtJ,...,LtJ-i} 



.(2j + l)7r ■(2k}n_ 



P {e " , e " 
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where Aq is the product of the non-zero eigenvalues of Dq. Let us define: 

gix, y)=P (e-^ e^^) , /(x, y) = log \g{x, y)\=^ \og{g{x, y)), 
then taking the logarithm of the ratio of the two determinants yields: 
I det Kl^ 



log 



detB 



logAo+ ^n^f 

j,fce{-LfJ,...,LfJ-i} 
(i,fc)7^{o,o) 



2j 2k\ , 1 



n ^ n J 



ee{-l,+l} 



(14) 



where 



A 



ee{-l,+l} 



The first term in ()14p . — logAo is a constant. The last sum is equivalent to 21ogn 
because, again by Theorem [8l (1, 1) is a zero of P of order exactly 2, and thus, for a 
certain c 7^ 0, 



/(-,0)=log 



1 + 



2 log n + bounded error terms. 



To finish the proof, we need to show that first sum is o(21ogn). We first use Taylor 
expansion to show that for {j,k) ^ (0,0), 



Aa.;(?i.?i).J,|!/(?i,?i)+o('i sup 

\n n J ax^ \ n n J xc\ '^''^ 2j+i 



dx"^ \ ^ n J 



(15) 



n ' n 



Let F be a small neighborhood of (0, 0) in [—1, 1]^ contained in a ball of radius 6 G (O, ^) , 
and define 



n 




n 


-2. 


, . . . , 


.2- 



1 



I ^2 



{j,k) :(^,^^)ey^ 



As long as (j, A;) is outside of Vn, the fourth derivative of / is uniformly bounded, and 
the discrete sum restricted to those indices {j, k) converges to the integral of the smooth 
function on the square [—1, 1]^ deprived from V. Thus this part of the sum stays 
uniformly bounded as n goes to infinity. 

We now need to control the contribution to the first sum of p4|) of the indices that are 
in Vn, i-e. relatively close to zero. Since (1, 1) is a double zero of P, we have for x and 
y close to zero the following expansions for g and its derivatives: 

gix, y) = P(e^-^ e'^y) = ax^ + 2(3xy + ^y" + 0(|| (x, y)f ), 



dx 



(x,2/) = 2(ax + /3y)+0(||(x,y)f), 



5^ 
dx^ 



(x,y) =2Q + 0(||(x,y)||), 
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with a, f5, 7 real numbers. Moreover, since (1, 1) is the only zero of P on the unit torus, 
the quadratic form 

q{x, y) = ax^ + 2(3xy + -fy"^ = {{ax + f3yf + (07 - , 

is definite. Without loss of generality, we can suppose that a is positive so that q 
defines a scalar product on M^. 

Let us now compute an expansion of the second derivative of / with respect to the first 
variable: 

52 f 9{x,y)Q{x,y) - (^{x,y)) 
^{x, y) = ^ 



2a' 



g{x,y)^ 

-{ax + Py)^ + {aj - (3'^)y^ 
{{ax + I3yf + (a7 - fi'^)y'^f 



:i + 0(||(x,y)||)), (16) 



if has the form ^j, then (fT6]) divided by 'n? becomes: 

-(aj + /?A;)2 + (a7-/32)fc2 



n 



2 dx"^ \ n 



a 



{{aj + (3kf + {a-f - (5^)k^y 



1 + 



n n 



(17) 

When we sum over (j, k) m Vn \ {(0,0)}, the contribution of the error term is finite, 
since for some constant C > 0, 



(j,fc)7^{0,0) 



-{aj + l3kf + (a7-/?2)A;2 



{{aj + l3kY + {a-i - l3'^)k'^Y 



X O 



i. h. 

n n 



1 



||(i,fc)||<<5n U ^ J 
{j,fc)^(0,0) 



which is bounded uniformly in n, as can be easily seen by estimating the sum using 
polar coordinates. 

Let us now take care of the main contribution. We see that the expression 

a^ ( -{aj + pkf + {aj - I3^)k'^ \ 
2 [({aj + Pk)^ + {aj - (3^)k^f ) ' 
is changed in its negative if we change (j, k) so that the following two quantities 

aj + pk and \/cry — P^k, 

are interchanged, i.e. if we replace {j,k) by its image by the symmetry with respect 
to the bisector of the two straight lines ax + /3y = and y'cry — JPy = along its 
orthogonal for the scalar product q. Unfortunately, this transformation does not in 
general send an element of to another couple of integers. But if we call {j',k') 
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the closest lattice point to the image of {j, k) by this symmetry, then the sum of the 
contributions of {j,k) and {j',k') is of order 0(||(j, A;)||~^) and thus is summable. If 
the neighborhood V we took around the origin is invariant under this symmetry, then 
the sum over all (j, k) in this neighborhood of ()17p is bounded as n goes to oo. 

The fourth derivative of / is a 0(||(x, y)||~^) and thus 



sup 

• 2j-l 22+1] 



dx* \ ' n J 



O 



which is summable over Z^. The contribution of the error term in (|15p is bounded 
as n goes to infinity. All the contributions from the first sum in the RHS of ()14p are 
bounded when n gets large. So we get finally that 



log 



detB 



< log 



detB 



det 



21ogn + 0(l), 



which ends the proof of Part 2 of Lemma [71 



□ 



5 Characteristic polynomial and the Laplacian 

Let G be an infinite Z^-periodic isoradial graph, and let Q be the corresponding Fisher 
graph. In this section, we define the critical Laplacian operator of the graph G, the 
Laplacian characteristic polynomial and the Laplacian spectral curve. We also recall the 
definition of the characteristic polynomial of the dimer model on Q. Part 1 of Theorem 
[8] states that the Laplacian spectral curve is a Harnack curve of genus with a unique 
point (1,1) on the unit torus, it is proved in Section [5Tl Part 2 of Theorem [8] states 
that the dimer spectral curve of Q and the Laplacian spectral curve of G are the same, 
it is proved in Section [5.21 

The critical weights for the Laplacian are tan(0) |Ken02| . The Laplacian matrix on G, 
corresponding to these weights, is defined by: 



tSLn{9uv) if M ~ f, 

- En'^u tan{0uu') iiu = v, 
otherwise. 



This naturally defines the Laplacian operator acting on functions on V{G). 

In a similar way to what we did for the Kasteleyn operator in Section r4.2.1l let A{z, w) 
be the Fourier transform of the Laplacian A. The Laplacian characteristic polynomial, 
denoted P/\{z,w), is the determinant of A{z,w). The set 

{{z,w) eC^ : Pa{z,w) =0} 
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is called the Laplacian spectral curve of G. 

Recall the definition of the characteristic polynomial P{z, w) of the critical dimer model 
on Q: ^ 

P{z,w) = det{k{z,w)), 

where K{z, w) is the Fourier transform of the Kasteleyn matrix K of the dimer model 
on Q. The main result of this section is the following: 

Theorem 8. 

1. The Laplacian spectral curve of G is a Harnack curve of genus 0. Its unique point 
on is (1, 1), and it has multiplicity 2. 

2. There exists a constant c ^ 0, such that 

P{z,w) = cPa{z,w). 

The proof of Theorem [8] is given in the next two sections. 

5.1 Laplacian spectral curve: proof of Part 1 of Theorem [H] 

In order to prove Part 1 of Theorem [HI we show in Proposition [10] that the Laplacian 
characteristic polynomial is equal, up to an explicit multiplicative constant, to the 
characteristic polynomial of the dimer model on the double of the graph G, denoted 
Gil, defined below. The graph Gd is isoradial and bipartite, so that by Kenyon and 
Okounkov [KO06j . we know that the dimer spectral curve of Gd is a Harnack curve of 
genus 0. 

Moreover by definition of the Laplacian, the sum of the columns of the matrix A(l, 1) 
is the zero vector, so that det(A(l, 1)) = -Pa(I) 1) = 0. Since a Harnack curve has or 
2 conjugate zeros on a torus of given radius, we know that (1, 1) is the unique double 
zero of Pa{z, w) on T^. 

The statement of Proposition [TU] is given at the end of this section, since it is more 
natural to prove it before stating it. 

5.1.1 The double of the graph G 

The double of the graph G, denoted by G^, is the bipartite planar graph constructed 
from G as follows: there is a black vertex in Gd for every vertex of G and every vertex 
of G* , and there is a white vertex in the middle of each edge of G (or each edge of G*). 
A white vertex w and a black vertex b are connected by an edge in Gd if the vertex of 
G or G* corresponding to b is adjacent to the edge corresponding to w, see Figure [H 
The double Gd is again a periodic isoradial graph, and its diamond graph is obtained 
from G^ by splitting each rhombus in 4 identical smaller rhombi. 



25 




Figure 6: A piece of G, together with the rhombi corresponding to the represented 
edges (left). The corresponding piece in the double Gd (right). The diamond graph of 
Gd is obtained by splitting every rhombus of G^ in four identical smaller rhombi. 



Note that by construction, the double of the dual graph G* is also Ga- 
in |Ken02] . Kenyon defines a bipartite complex Kasteleyn operator /C on Gd as an 
infinite matrix whose rows (resp. columns) are indexed by white (resp. black) vertices 
of Gd- If w is a white vertex and b a black vertex of Gd, then ICw,b = when they are 
not neighbors. When they are adjacent to each other, then 

/C^,b = : = 2sme^be' 2 , (18) 

where and e*^ are the two unit vectors based at w representing the two sides of the 
rhombus containing the edge wb, see Figure [7] below. 




b 



Figure 7: A rhombus of the diamond graph of Gd- It contains an edge between a white 
vertex w and a black vertex b- The sides of the rhombus are represented by the two 
unit complex numbers e*" and e*'^. 



5.1.2 Kasteleyn operator on Gd, incidence matrix and the Laplacian 

Consider any orientation of the edges of G- Then, for every white vertex w of Gd 
(corresponding to an edge of G), and for every vertex b of G, the coefficient of 
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the oriented weighted incidence matrix M*^, is defined by: 

{0 if tfie edges corresponding to w is not adjacent to b 

— ^/tanO if the edge corresponding to w is oriented away from b 
Vtan 9 if the edge corresponding to w is oriented towards b, 

where 9 is the half-angle of the rhombus of containing the edge of G corresponding 
to w. M*^*, the incidence matrix of G* is defined as above, with the orientation of the 
edges of G* induced by that of G. 

In order to avoid confusion, let us denote by A*^ and A*^* the Laplacian on G and G* 
respectively. From the definition, we then have: 

{M^)'^M^ = ISp, [M^')'^M^* = ISP\ {M^)'^M^* = 0. (19) 



Following Kenyon |Ken02j . we now relate the Kasteleyn matrix /C, the incidence ma- 
trices and M^* and the Laplacians A*^ and A*^*. Define a diagonal operator A 
on white vertices of G^ by: 

2vsm 9 cos 9 

where e^^ is the unit vector pointing in the same direction as the oriented edge of G 
corresponding to w, and 9 is the half-angle of the rhombus containing this edge. If we 
order the columns of /C by listing first the black vertices corresponding to vertices of 
G, and then vertices of G*, the product AKL has the following form: 

AK. = (M^ iM^') . 

Using (fT9]l . this yields: 

K^A^AK = (AIC)HAK) = (M° 'M°l = (f ^o-) ■ (21) 

All the above operators commute with translations of 1?. Equation (j21|) induces a 
relation on the Fourier transforms, which are finite matrices with polynomial entries in 
the variables G T'^: 

BItl£ = idA^AlC = 

Taking the determinant on both sides, and using the fact that: 

1 




del A 



,,/^ N r.,^ N 2A/sin 9e cos 9^ ' 

u'6y(Gdi) ee£;(Gi) 
w white 
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yields the following relation: for every {z,w) £ T^, 



n 



1 



4 sin 9p cos Op 



det K:{z,w) = detAG{z,w)detA'^*{z,w). (22) 



The next step consists in relating det A*^(z, tt;) and det A'^* (z , w) . This is done by 
giving a combinatorial interpretation of det A'-' {z , w) = Pj\{z,'w) in terms of cycle- 
rooted spanning forests on Gi. This is very similar to the link between the determinant 
of a principal minor of the Laplacian of a graph and spanning trees on this graph, given 
by Kirchhoff's matrix-tree theorem |Kir47j . 



5.1.3 Combinatorial interpretation of Pa(z,w) and conclusion of Part 1 

A cycle-rooted tree T of Gi is a connected subgraph of Gi with a unique cycle, which 
must be non trivial, i.e. not contractible to a point on the torus on which Gi is drawn. 
The homology class of T is that of its non trivial cycle in Z^, which is well defined 
up to a sign since the cycle is not oriented. A cycle-rooted spanning forest of Gi is a 
collection of disjoint cycle-rooted trees covering every vertex of Gi. Let us denote by 
!F the set of cycle-rooted spanning forests of Gi . 

Lemma 9. The Laplacian characteristic polynomial is the following combinatorial sum: 
Pa(^,'"^) = ^ Jl [W tan 06 \/z^wy - ^Jz-'^w-y ^ | , (23) 

where ±{x,y) £ 1? is the homology class of the cycle-rooted tree T. 

Proof. We follow the same strategy as for the proof of the matrix-tree theorem (see e.g. 
|Big93| ). Since we don't use the dual graph G* in this argument, in order to improve 
readability, we drop the superscript G in the operators A*^ and M*^, and their Fourier 
transform. We start from the Fourier transform of the first relation in (I19p 

A = Wm. 

Recall that the rows of M are indexed by edges of Gi , and its columns are indexed by 
vertices of Gi. Let us use Cauchy-Binet's Theorem to express the determinant of A as 

^ I 2 

det A(z, w) = det (Aft ( z,w)^)det{M{z,w)s) = ^ det(M(2:, w)^) . 

S^E{Gi) S<zE{Gi) 
\S\=\V{G^)\ |5| = |V(Gi)| 

The matrix M{z,w)s is the incidence matrix of the graph Gs obtained from Gi by 
taking only the edges in S. If Gg contains a trivial cycle, then we can find a function 
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with support on this cycle in the kernel of Ms- Therefore the only contributions come 
from graphs, whose connected components have a number of edges equal to the number 
of vertices, and each have a non trivial cycle, i.e. a cycle-rooted tree. All these non 
trivial cycle are parallel, i.e. they have the same homology, up to a sign. 

After a possible reordering of the rows and of the columns of the matrix M{z, w)s, we 
can assume that it is block diagonal, each diagonal block corresponding to a connected 
component. The square root of the weight of each edge in S can be factored out from 
the determinant from each row of M{z,w)s. The determinant of each block can be 
evaluated recursively by expanding it along columns corresponding to vertices of degree 
1, i.e. leaves of the tree, collecting a product of terms of modulus 1 along the branches. 
The evaluation of the determinant reduced to the cycle gives two terms of modulus 1, 
and the ratio between these two terms is —z^w^. □ 



A direct consequence of this lemma is: 



detAS{z,w)=^ Yl tan^e detA^*(z,w). (24) 

Indeed, the dual T* of a cycle rooted spanning forest T on Gi, consisting of edges e* 
such that e T, is a cycle-rooted spanning forest on G^. Moreover, since the weight of 
an edge for the Laplacian is the reciprocal of its dual edge, 

tan^e* = tan ^— — 9^ = (tan^g)"^, 

the weight of T* is equal to the weight of T divided by IleG-BCGi) tan^e- 
Combining ^ and dM} yields: 



Proposition 10. There exists a complex number ( of modulus 1, such that for every 
{z,w) G €?, 

det AG(z,'u;) = C JJ 2 cos detK,{z,w). 
Vee£;(Gi) / 

This ends the proof of Part 1 of Theorem [8l In the next section we prove Part 2, namely 
that the dimer characteristic polynomial P of ^ is equal to up to a multiplicative 
constant. 



5.2 Dimer and Laplacian characteristic polynomial: Proof of Part 2 
of Theorem [H] 

In order to prove part 2 of Theorem [HI we first show in Lemma [TT] that P{z, w) is 
a divisor of P^{z,w), and then in Lemma [T2] that the Newton polygon of P{z,w) is 
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included in the Newton polygon of PAizjw). Combining these two results yields that 
there exists a constant c 7^ 0, such that 

\/{z,w)gC'^, det K {z , w) = c det A{z , w) . 

Lemma 11. P{z,w) is a divisor of Pa{z,w). 

Proof. Let T^'^{G) be the space of (z, ii;)-quasiperiodic functions 

T^'^iG) = {/ : Z2 ^ C^(^) : V(x,y) e Z\ f{x,y) = z'^w'^ f (0,0)} . 

A function of this space is entirely determined by its value on a fundamental domain. 
This space is thus naturally in bijection with C^^*^^^ by the following projection: 

/^/(0,0) 

Since A is translation invariant, we have for any function / G J^^''^{G) that 

V(x,y)eZ2, ^veViGi), {A^ f)ix,y), = z-^w-y (AG{z,w)f (0,0)) ^. 

A'-'{z,w) can thus be interpreted as the matrix of the restriction of A*^ to J-^''^{G) in 
its canonical basis. 

By |Ken02l IMerOlj . the discrete exponential function on G with parameter A G C is 
harmonic [i.e in the kernel of A) and (2;(A), t(;(A))-periodic for 

-(A)= n n (25) 

e=exp{ia)e£{G«) e=exp(ja)e£;(G«) 

The projection of this exponential function on Gi is a vector in the kernel of the matrix 
A^{z{X),w{X)). Thus {z{\),w{X)) is a point of the spectral curve. The Laplacian 
spectral curve is Harnack, and thus irreducible. Therefore A 1— > {z{X),w{X)) gives a 
parametrization of the whole curve, except maybe some isolated points. 

But for any A G C, we can also construct a (z(A), ti;(A))-periodic function f\ on vertices 
of Q such that 

Vr; G V{G), Kfx{v) = 0. (26) 

See |BdT09] for the details of the construction. Therefore, the matrix K{z{X),w{X)) is 
not invertible, and {z{X),w{X)) is a point in the spectral curve {{z, w) £ C'^ : P{z, w) = 
0}. This shows that Pa is a divisor of the characteristic polynomial P. □ 



30 



We now have to show that the quotient of these two polynomials is a constant. The 
idea is to compare the degrees. Recall that the Newton polygon N(P) of a polynomial 
P{z, w) is defined as the convex hull in M2 of 

{{x, y) G : z^w^ is a monomial in P{z, w)^ . 

We show that the following inclusion holds: 

Lemma 12. 

N{P) C A^(Pa). 
A key ingredient to the proof of Lemma [12] is the following: 

Lemma 13. For every monomial of P{z,w) = det(^K{z,w)^ , one can construct a 
configuration of non intersecting, non trivial oriented loops on the toroidal isoradial 
graph Gi, with homology equal to the degree of the monomial. 

Before proving Lemma [T3l let us end the proof of Lemma [T2l 

Proof of Lemma [TB. Take an extremal point (xq, yo) of the Newton polygon of P{z, w), 
then z^^w^" is a monomial of P{z,w). Construct a family of n non trivial loops ac- 
cording to Lemma [T31 Then each loop has homology ±{x,y) = ib(^,^), where 
d = gcd(xo,yo)- From this set of loops, one can grow a cycle-rooted spanning forest F 
by gluing new edges step by step. This forest has n trees, and the weight of the forest, 
computed as in (I23p is a polynomial in z and w, with monomials of degree in 

DF = {-n,...,n}ix,y). (27) 

The set of possible degrees Dp in the weight of a spanning forest constructed from this 
family of loops contains (xo,yo)- Note that the terms with extremal degree come with 
positive coefficients. It follows that the whole set Dp is included in the Newton polygon 
of N{P/\). It might happen that the monomials corresponding to the extremal points 
oi Dp are canceled "accidentally" by the contributions of other forests F' , but the cor- 
responding monomials in these contributions have to come with a negative coefficient, 
implying that the extremal degree of these terms is larger, and the corresponding set 
of possible degrees Dp' contains Dp. Thus the extremal point (xcyo)) ^-s an element 
of Dp, belongs to iV(PA)- 

Since N{P/^) contains all the extremal points of N{P), it contains all the Newton 
polygon N{P) by convexity. □ 

Proof of Lemma In order to construct the family of loops on Gi , we need some 
combinatorial description of det[K{z, w)) . The matrix K{z, w) is not skew-symmetric, 
but only skew-Hermitian, because of the complex entries z and w. Thus, the deter- 
minant of K{z, w) is not the square of a combinatorial sum over dimer configurations. 
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that would correspond to the superposition of two independent dimer configurations. 
However, we can stih expand the determinant as a sum over permutations, and group 
terms according to the support of their cycles. 



det{K{z,w))= sgn{a) [] {K{z,w)) ^^^^ = 

(r^) partition i 7 cycle veT^ 

of y (6-1) supp7=r^ 
One can make a few remarks about this expression: 

1. Since K{z,w) is an adjacency matrix of Qi, the only non zero contributions to 
this sum will come from partitions of Qi into disjoint circuits or doubled edges. 

2. If is a circuit, the sum over cycles with support contains two terms corre- 
sponding to the natural cycle 7 sending a vertex to its successor on the circuit, 
and its inverse 7"^. When is a doubled edge, the sum contains a single term, 
representing the transposition of the two extremities, which is equal to its inverse. 

3. The degree of {z, w) of the product n?;er(; {^i^^ ^'))v'yv homology in 

of the oriented cycle 7 around the torus on which ^1 is drawn. Two cycles 7 and 
7~^ have opposite homology. 

4. Because K{z, w) is skcw-Hermitian, the coefficients of the monomials for 7 and 
7~^ will be equal if |r^| is even, and are negative of each other if |r^| is odd. 

5. It follows from the two previous remarks that trivial loops of odd length do not 
contribute (the terms from 7 and 7"^ cancel each other). Recall that on Qi there 
are triangles around each decoration, which are trivial cycles of odd length. We 
will for the moment only cancel out the contributions of configurations containing 
a cycle around a small triangle of a decoration. 

6. The weight corresponding to a transposition is a positive number. 

7. Since they cannot cross, cycles of non trivial homology are parallel: their homol- 
ogy class in are equal or opposite to one another. 

We say that a cycle on ^1 visits a decoration if the sequence of its oriented edges 
contains successively: 

• a long edge {i.e. an edge coming from Gi) oriented toward the decoration, 

• a sequence of distinct edges of the decoration. 
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Figure 8: A configuration of two cycles visiting a decoration in different directions: in 
the initial state (left) and after rerouting (right). 



• another long edge attached to the decoration, but oriented outward the decora- 



A cycle can visit a decoration in two possible directions (clockwise or counterclockwise), 
depending on which sector of the decoration is taken to connect the two long edges. A 
cycle can visit several times the same decoration, but always in the same direction. 

Take a monomial a^^yZ^w^ of P{z^ w) = det K{z, w). Consider a configuration of cycles 
on Qi contributing to this monomial in the expansion of the determinant (j28|) . 

Consider the set of cycles visiting a particular decoration. If these cycles do not all 
visit the decoration in the same direction, we can construct just by elementary moves 
around the decoration a cycle configuration also contributing to the same monomial, 
but with a weight equal to the exact negative of the weight of the initial configuration. 
Pairing these configurations gives thus a zero contribution to (j28p . 

The construction is done as follows: since not all the cycles visit the decoration in 
the same direction, then one can find two neighboring cycles with opposite direction. 
The sector of the decoration separating the two loops is filled with doubled edges with 
weight 1, since we already threw away configurations with cycles around triangles. Now 
reroute the cycles around the triangles, and switch all the doubled edges (see Figure [8]). 
The two cycles either increased or decreased their length by 1 unit, preserving the 
global signature. The condition on the Kasteleyn orientation implies that the product 



where ui,U2,U3 are vertices of the same small triangle. Consequently, when the cycles 
are rerouted, the product of the entries of the matrix K for one of the two cycles will 
stay the same, whereas the one for the other cycle will be changed in its negative. 
Therefore the total weight of this new cycle configuration is exactly the negative of the 
weight of the initial configuration we started with. 



tion. 
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Having paired configurations of cycles witli anniliilating weiglits, we are left with con- 
figurations where every decoration is visited by cycles in the same direction. If there is 
a decoration visited by two non trivial cycles (there cannot be more, for obvious topo- 
logical reasons), then they must be parallel. If they have the same homology class in 
Z^, they visit the decoration with different orientations. Therefore, they have opposite 
homology, and thus they do not contribute to the homology class of the cycle configu- 
ration. We can erase these two non trivial loops, and search for other decorations on 
which we can reiterate this operation. 

Once that around every decoration, there is at most one non trivial cycle, one can 
throw away all the trivial ones as well as the double edges, and retract every decoration 
to a single vertex. This would lead to a configurations of non intersecting loops on Gi 
with homology {x,y), which is exactly what we wanted to construct. 

□ 
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